A Lower Bound on Neutrino Mass and Its Implication on The Z-burst Scenario 



Kwang-Chang LaQ 
Department of Physics, National Taiwan University, Taipei, Taiwan. 

Pisin Cheilll 

Kavli Institute for Particle Astrophysics and Cosmology, 
Stanford Linear Accelerator Center, Stanford University, Stanford, CA 94309 
(Dated: February 5, 2008) 

We show that the cascade Umit on ultra high energy cosmic neutrino (UHECi^) flux imposes a 
lower bound on the neutrino mass provided that super-GZK events of ultra high energy cosmic rays 
(UHECRs) are produced from Z-bursts. Based on the data from HiRes and AGASA, the obtained 
neutrino mass lower bound violates its existing cosmological upper bound. We conclude that the Z- 
burst cannot be the dominant source for the observed super-GZK UHECR events. This is consistent 
with the recent ANITA-lite data. 

PACS numbers: 



I. INTRODUCTION 

Big bang cosmology predicts the existence of both 
cosmic microwave background (CMB) and cosmic 
neutrino background (C:^B). Ultra high energy cosmic 
protons are expected to interact effectively with the 
CMB photons, predominantly through the photopion 
production at A-hyperon resonance, and would lose their 
energies rapidly with the attenuation length around 
50Mpc. As such the ultra high energy cosmic ray 
(UHECR) spectrum is predicted to exhibit a cutoff — 
the so called GZK cutoff [3, 01 — around the threshold 
energy 4 x lO^^eV. While observations from the HiRes 
experiment is consistent with the notion of GZK cutoff 
the AGASA data appears to suggest the opposite 
This leads to many speculations as to whether the 
GZK cutoff really exists, and if not what is the nature 
of these super-GZK events. 

Existing models for super-GZK UHECRs are usually 
categorized into top-down and bottom-up scenarios. 
The top-down scenario assumes the existence of super- 
massive exotic elementary particles based on theories 
beyond the standard model. The major challenge 
of this scenario lies in the demand for a fine-tuned 
decay and/or annihilation rate and the lack of physical 
evidence for their existence. On the other hand, the 
bottom-up scenario, which assumes ordinary particles 
as the UHECRs, faces the challenge of providing an 
effective mechanism to accelerate particles to ultra high 
energies. Even if an effective "cosmic accelerator" can 
be identified, the issue of cosmic transport dictated by 
the GZK mechanism remains, as there does not seem to 
exist identifiable sources within our local super cluster 
50Mpc) for the detected events. 



To circumvent this difficulty, it was suggested that the 
Z-burst, the resonant annihilation of the ultra high en- 
ergy cosmic neutrino (UHECz^) with the Ci^B into a Z 
boson and its subsequent decay into ultra high energy 
protons that occurs within our local super cluster 

can account for UHECRs beyond the GZK-cutoff |i[l3. 
With its mean-free-path comparable to the present Hub- 
ble radius, the UHEC;^ serves as a cosmic messenger that 
can avoid the GZK proton attenuation problem without 
invoking particle theory beyond the standard model. For 
the Z-burst to happen, the UHECi^ must be at a resonant 
energy. 
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which depends on the neutrino rest mass, mj/. Here Mz 
denotes the mass of the Z boson. If the Z-burst mech- 
anism is indeed responsible for the observed UHECR 
super-GZK spectrum, then there must exist a constraint 
on the neutrino mass via the above relation. 

Based on the Z-burst scenario, two groups have 
derived bounds on neutrino masses using AGASA data 
with different strategies. Fodor, Katz and Ringwald 
deduce the Z-burst spectrum from the AGASA data 
by parameterizing the transition from the non-burst to 
the burst component near and above the "ankle" of 
the UHECR spectrum. Gelmini, Varieschi and Weiler 
derive their bound by requiring the Z-burst not to 
overproduce non-observation events beyond the AGASA 
end-point energy. Our strategy, instead, is to invoke an 
upper limit on the UHECj/ flux so as to obtain an upper 
bound on the required resonant energy Eres, which can 
in turn be translated into a lower bound on the neutrino 
mass. 
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In this paper we derive a lower bound on the neu- 
trino mass based on the assumption that Z-burst mech- 
anism saturates the observed UHECR super-GZK spec- 
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trum. Our deduced neutrino mass lower bound, however, 
turns out to be higher than the existing upper bound 
deduced from cosmological considerations. We thus con- 
chidc that the Z-burst mechanism cannot be responsible 
for the supcr-GZK UHECR spectrum. Our conclusion 
agrees with that from the recent ANITA-lite experiment 



that for a single neutrino flavor. We further assume that 
UHECz^ fluxes are the same for neutrinos and antineutri- 
nos. By definition. 



dr / dEF{E,r) 
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II. UHECR FLUX AND Z-BURST 

Assume that all observed super-GZK UHECR proton 
events are induced from Z-bursts. The observed super- 
GZK proton flux must be smaller than the total Z-burst 
proton yield in the universe since there must be events oc- 
curred outside our local GZK-sphere 50 Mpc) which 
could not reach the Earth. Furthermore, in order for 
Z-burst events to saturate the observed super-GZK spec- 
trum, it is inevitable that they are oversupplied since 
there must be some Z-burst protons that are generated 
at energies below the GZK-cutoff. Therefore, 



Tobs ^ jZ 

^p,>GZK ^ ^p\u ■ 



(2) 



where Ip^^QZK total observed super-GZK proton 

flux with energy exceeding the GZK-cutoff and I^^^ 
the total proton flux from Z-bursts, both in units of 



Here F{E, r) is the UHEC;/ flux at energy E and 
distance r from the Earth, ni,{r) is the number density 
of the CvB, (T,yf^{s) the neutrino-antineutrino cross 
section at s = 2myE^ Br(Z hadrons) the branching 
ratio, and S,p+n\u the multiplicity of nucleons per Z-burst. 

For completeness, our integration should include all 
neutrinos and Z-burst events in the universe. Such a 
treatment tends to be over-conservative as the protons 
deduced from Z-bursts outside of our local GZK-sphere 
may hardly survive. The complete but ultra-conservative 
treatment is discussed in the appendix. A physically rea- 
sonable yet much simplified calculation can be carried out 
by neglecting the contributions outside of our local GZK- 
sphere. This amounts to replacing the maximum distance 
Rmax in our integration by the radius of our local GZK- 
sphere {Rgzk ~ 50Mpc). As the distance under consid- 
eration is much more local, all the r-dependence can be 
ignored: 



Though observations cannot completely rule out 
the possible contribution to the super-GZK UHECR 
spectrum by UHE photons, experimental data [T^ 
suggests that protons saturate the super-GZK flux, 
^P^>GZK ^ ^>GZK^ ^t 2(7 confidence level. Then in terms 
of the total observed UHECR flux, Eq.® can be written 
as 



johs ^ jZ 
^>GZK ^ ^p\i, ■ 



(3) 



The AGASA experiment has accumulated 57 events 
above 4 x lO^^eV with a total exposure of ~ 4 x 
10^'^cm^ ssr 15]. This translates into an observational 
super-GZK flux. 
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1.43 X 10~^^cm~2g-igj.-i_ 



(4) 



It can be shown that a fitting spectrum with a power-law 
index —2.78 J^J reproduces the above flux. 



III. Z-BURST YIELD 

Now we deduce the total Z-burst proton yield within a 
relevant cosmic volume. Solar and atmospheric data on 
neutrino oscillations indicate that the oscillation lengths 
are much shorter than the solar distance. So for cosmic 
neutrinos their population among the 3 flavors should be 
equalized. The total UHEC;^ flux is thus simply 3 times 



Ip\i' = 35p+„|^i?GZi^f^l/Br(Z -> hadrons) 



poo 
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where n,j{r) = rii, — 112cm ^ is the neutrino- 
antineutrino number density per flavor at present. 

The UHECi' flux is commonly assumed to follow a 
power-law energy spectrum 



F{E) = FoE-", 



(7) 



with Fq being the normalization factor. 

Using E = sEres/M'^, we can write the energy inte- 
gration as 

/ dEF{E)a„^{s ^2mE) 
Jo 

/■°° ds 

= Eres / — jF(s£;,es/Af|)a,p(s). (8) 
Jo 

As the neutrino-antineutrino annihilation cross section 
is sharply peaked at the Z-resonance, it acts essentially 
like a J-function in the integration over the energy of 
the UHECi/. We therefore introduce the energy-averaged 
cross section 0, ^| 

/ds 
— 2 fT,p(s) = 2^V2Gi. = 40.4nb, (9) 
z 
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which is the effective cross section for all neutrinos within 
the resonance range {Eres{'^—^ z /Mz), Eres{^+^ z /Mz)) 
and simplify the integration of Eq. ^ as 

f°° ds 

JQ ^"Z 



Putting everything together we find 

M2 



(10) 



RGZKn„^p+„\^FQ \^^— j 

X (cr^p)Br(Z ^ hadrons). (11) 

The experimental data ^3 gives the branching ratio 

Bt{Z hadrons) = (69.89 ± 0.07)%. (12) 

The final proton multiplicity per Z-burst was calculated 
by Fodor, Katz and Ringwald(see |1^ and references 
therein) as a function of the proton momentum distri- 
bution and by Gelmini, Varieschi and Weiler jl^ us- 
ing the event generator PYTHIA They obtain 
^p^n\i^ = 2.04 and 1.6, respectively. We take the former 
value in this paper. 



IV. CASCADE LIMIT AND NEUTRINO MASS 
BOUND 

We now invoke the cascade limit to constrain the 
VSECv flux 0, I23. This is permissible due to the 
fact that neutrino productions must always be accompa- 
nied by photons and electrons. The cascades are induced 
while these photons or electrons interact with low energy 
background radiations such as the CMB in extra galactic 
space and the infrared radiation inside the galaxy. The 
photons so induced would further cascade and eventu- 
ally pile up in the energy range of lOMeV-lOOGeV with 
a spectrum oc E~'^, which is consistent with the EGRET 
observation The estimated average energy density 

in this range is ^'-'egret ~ 2x 10~®eV/cm'^. This provides 
an upper bound on the UHE neutrino flux, 



E'^F{E) < —LJcas < -j— WEGRET 



(13) 



To be prudent, we do not assume the exact value of a = 2 
for the power-law index, but instead leave a as a free 
parameter, knowing that its value should be close to 2. 
Thus the parameter Fq can be substituted with an upper 
bound as follows: 



Fq < -^luegket / E:^ 
Fq < j^ujegket/E:^ 



2-a 
rain 
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a > 2, 
a < 2, 



(14) 



where E^ax a-nd Ej^in £^re the maximum and minimum 



energies of the UHECi^ spectrum. 



Implementing the cascade limit condition, and insert- 
ing all the relevant physical quantities discussed in the 
previous section, Eq.Q becomes 



E^, 



> 28. 7(-^)"-^ eV , a>2, 

> 28.7(:|^)2-"eV , a<2. 



(15) 



Note that in this expression not all the nii, dependence 
were grouped to the LHS, as E^-es clearly depends on 
•niij. Nevertheless this expression has an advantage in 
that EreslEmin > 1 and Emax/Eres > 1 by definition. 

An explicit lower bound can be obtained by moving 
all the nil, dependence to the LHS. We then find 



> \{A^^eI-^^MI) , a>2, 

1 1 2^ 

> -{A~E,l-alMl) , a<2, 



(16) 
(17) 



where 
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(18) 



The mass bound is dependent on the power-law in- 
dex a and the values of Emax or Emin- Our limited 
knowledge on the UHECz^ renders large uncertainty in 
the determination of Emax and Emm- One thing which 
is certain, however, is that the resonant energy must lie 
in between Emax and Emin for the Z-burst to happen. 
Ea. (|15|) indicates that the minimum value of our bound 
corresponds to the situation where either Emax or Emin 
equals Eres, or a = 2. Since we should look for the low- 
est possible lower bound, we put Eres = Emax = Emin 

in our estimate and arrive at our neutrino mass lower 
bound 

m,>28.7tllieY {Rmax ^ Rgzk 50Mpc), (19) 
where the error comes from fitting the AGASA data . 

Recent WMAP ^ measurement of the CMB fluc- 
tuations has deduced a strong upper limit on neutrino 
masses, Tjim^. < 0.69eV. Since any single neutrino mass 
~ 0.04eV implies a near mass-degeneracy for all three ac- 
tive neutrinos, one concludes < 0.23eV. Two analy- 
ses mm which include data from WMAP, 2dF, SDSS, 
and galaxy cluster surveys have arrived at a bound of 
^ 0.7eV. Another analysis using CMB and 
LSS data gives EiTOi^. < leV, but finds a stronger bound 
^if^ui ^ 0.6eV when priors from supernova data and 
Hubble Key Project are included. These newer results 
are close to the original WMAP bound. All these anal- 
yses converge to a cosmological upper bound of < 
0.23eV, which is 2 order of magnitudes smaller than the 
lowest possible lower bound we have derived. 
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V. IMPLICATION 

Our derivation is based on two assumptions: the 
saturation of the observed super-GZK UHECR flux by 
the Z-burst mechanism and the cascade upper Umit on 
the maximum UHEi/ flux. 

Since the cascade hmit is deduced from the cascades 
of photons accompanying the neutrino production, it is 
vahd for aU sources. Not only astrophysical accelerators 
(e.g. GRBs, AGNs, SNs, etc.) but also top-down 
sources, such as topological defects, superheavy X 
particles, dark matter, etc., are all contributing to this 
limit as long as the photons are co-produced along-side 
with neutrinos. It is generally believed^^j ^^^.t the 
cascade limit on UHEi/ flux is quite robust. 



Even under the most conservative assumption, 
Emin = E„iax = E^es, i-G- that the UHEC;/ spectrum is 
a delta function, our deduced neutrino mass lower bound 
is more than 120 times larger than the existing upper 
bound. We therefore conclude that the Z-burst scenario 
cannot account for the observed super-GZK UHECR 
flux. Assuming all the parameters are fixed, our neutrino 
mass lower bound can be lowered if we allow Z-burst to 
contribute only partially to the observed UHECR flux. 
Based on our values we can conclude that the Z-burst 
contribution to UHECR cannot be more than ~ 1% 
within our local GZK-sphere. The recent ANITA-lite 
experiment (13.j indicates that the Z-burst can at best 
contribute ~ 10% to the UHECR spectrum, which is 
consistent with our conclusion. 

Analogous to the GZK process, Z-burst is one of the 
few robust cosmic interaction processes based on the 
standard model of particle physics. With a mean-free- 
path comparable to the present Hubble radius, it pro- 
vides much hope to resolve the existing challenge of the 
bottom-up scenario. Our negative conclusion on its vi- 
ability as a solution to the super-GZK puzzle seems to 
force us back to the original dilemma. If the excess super- 
GZK flux is found to be real, the need for a solution 
remains. 



where ilm ~ 0.3 and ~ 0.7 are the present matter 
and dark energy densities in units of the critical density, 
respectively. 

According to the big bang cosmology, the number den- 
sity of Cz^B is n^(z) — n^{l + z)'^cm^^, where rii, — 
112cm~'^ is the neutrino-antineutrino number density at 
present. The UHEC;^ flux is now assumed to follow a 
power-law energy spectrum with cosmological evolution 
of the source included [1, and can be parameter- 

ized as 

F{E, r) = F{E)f{z) = FoE-^fo{l + zf , (A2) 

where /g is the normalization factor for the source evo- 
lution determined by the condition 



f{z)dz = 1 , 



(A3) 



with Zmax being the redshift of the most distant source. 

A more sophisticated distribution function has been in- 
troduced, based on the star formation rate and the GRB 
site distribution 31] 



/(2) = /o 



l + ai 



(1 + z)-"2 + ai{l + zY'i 



(A4) 



where ai,a2 and 03 are fitting parameters. It can be 
shown that our neutrino mass bound is insensitive 
between these two choices of distribution functions and 
for simplicity we will invoke Eq.(8) in our subsequent 
discussion. The separation of variables allows us to 
carry out the UHECi^ energy and distance integrations 
independently. 

The integration over propagation distance involves sev- 
eral sources of cosmic evolution. 



cm 



n,{l+zff{z)dz 
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APPENDIX A: A COMPLETE BUT 
ULTRA-CONSERVATIVE VERSION 

For completeness, the integration over distance must 
be carried out to include the cosmological evolution of 
the CvB number density and the UHECi^ flux. Since the 
dependence on the propagation distance r of the UHEC:/ 
fluxes and the C^'B number density can be expressed 
more straight- forwardly in terms of the redshift parame- 
ter z, we make the following change of variables 



dr 



(1 + z)^/nm{i + z)^ + nA 



dz, (Al) 



where simple analytic result is not readily attainable. Ta- 
ble U displays values of f{zmax,P) for selected choices of 
Zmax and p. It is clear that f{zmaxj/3) is reasonably in- 
sensitive to Zmax and /?. 



TABLE I: Values of f{zmax,P) for selected Zmax and j3. 





-3 


-2 


-1 





1 


2 


3 


1 


1.47 


1.54 


1.62 


1.69 


1.76 


1.82 


1.88 


Zmax 2 


1.65 


1.81 


2.00 


2.18 


2.35 


2.49 


2.59 


3 


1.73 


1.98 


2.27 


2.56 


2.81 


2.99 


3.12 
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With the knowledge of f{zmzx,l3), the Z-burst yield, 
Eq.(|SJl, becomes 



Ml \ 



Ip\u = cH„'^n^f{zjnax,P)ip+,i\,,Fo(—-^) 



X ((Ti,p)Br(Z hadrons) . 



(A6) 



Correspondingly, the lower bound on neutrino mass is 
changed to 



> 1.25/(z„„„/3)-i(^)"-2 



eV 



a > 2, 



m„ > 1.2b f{Zmax,fi) ^( 



-i^^)2-"eV, a<2.(A7) 



while Eas. p6|l and remain the same form with 



A = L 
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Cp+n|^(c^^p)Br(Z hadrons) 
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The mass bound is now dependent on the numerical 
values of the evolution factor f{z„iax,P) as well (See 
Table A gain we look for the lowest possible lower 
bound by putting E^es = Emax = Emin in our esti- 
mate. As we have shown, f{zmax,P) is of the order 1 
and is insensitive to z^ax and (3. With the choice of 
fizmaxjP) = /(3, 0) = 2.56, we arrive at our neutrino 
mass lower bound 



m.>0.49;™j4eV 



(i?„ 



cHo'). (A9) 



The former error comes again from the AGASA data 
and the latter from the uncertainty of the evolution 

factor f{Zmax,P)- 



TABLE II: Values of mass lower bound for selected Zmax and 
l3 with a — 2 energy spectrum. 
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-3 
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3 


3 


Zmax 


1 


2 


3 


1 


2 


3 


1 


2 


3 


f (^Zmax , /5) 


1.47 


1.65 


1.73 


1.69 


2.18 


2.56 


1.88 


2.59 


3.12 


mass bound 


0.85 


0.76 


0.73 


0.74 


0.58 


0.49 


0.68 


0.48 


0.40 
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